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ABSTRACT 

Virial black-hole mass estimates are presented for 12698 quasars in the redshift interval 
0.1 < z < 2.1, based on modelling of spectra from the Sloan Digital Sky Survey (SDSS) 
first data release . The black-hole masses of the SDSS quasars are found to lie between 
~ 1O 7 M and an upper limit of ~ 3 x 1O 9 M0, entirely consistent with the largest 
black-hole masses found to date in the local Universe. The estimated Eddington ratios 
of the broad-line quasars (FWHM > 2000 km s _1 ) show a clear upper boundary at 
Lboi/LEdd — lj suggesting that the Eddington luminosity is still a relevant physical 
limit to the accretion rate of luminous broad-line quasars at z < 2. By combining 
the black-hole mass distribution of the SDSS quasars with the 2dF quasar luminosity 
function, the number density of active black holes at z ~ 2 is estimated as a function 
of mass. In addition, we independently estimate the local black-hole mass function 
for early-type galaxies using the M^h — ° an d — Lb u i ge correlations. Based on 
the SDSS velocity dispersion function and the 2MASS if —band luminosity function, 
both estimates are found to be consistent at the high-mass end (M^ > 1O 8 M0). By 
comparing the estimated number density of active black holes at z ~ 2 with the local 
mass density of dormant black holes, we set lower limits on the quasar lifetimes and 
find that the majority of black holes with mass > 10 8 5 M are in place by ~ 2. 

Key words: black hole physics - galaxies: active - galaxiesmuclei - quasars: general 



1 INTRODUCTION 

The discovery that supermassive black holes are ubiquitous 
among massive galaxies in the local Universe indicates that 
the majority of galaxies have passed through an active phase 
during their evolutionary history. Moreover, the strong cor- 
relation observed between black-hole and bulge mass (Kor- 
mendy & Richstone 1995; Magorrian et al. 1998; Gebhardt 
et al. 2000a; Ferrarese & Merritt 2000) indicates that the 
evolution of the central black hole and its host galaxy are 
intimately related. Consequently, it is clear that studying 
the evolution of quasar black-hole masses will provide cru- 
cial information concerning the evolution of both quasars 
and massive early- type galaxies. 

Within this context the last few years have seen re- 
newed interest in the possibilities of estimating the central 
black-hole masses of active galactic nuclei (AGN). The ma- 
jor impetus for this has been the results of the recent rever- 
beration mapping programmes carried out on low-redshift 
quasars and Seyfert galaxies (Wandel, Peterson & Malkan 
1999; Kaspi et al. 2000). The measurements of the broad-line 
region (BLR) radius produced by these long-term monitor- 
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ing programmes have allowed so-called virial black-hole mass 
estimates to be made for 34 low-redshift AGN (Kaspi et al. 
2000). The principal assumption underlying the virial mass 
estimate is simply that the dynamics of the BLR are domi- 
nated by the gravity of the central supermassive black hole. 
Under this assumption an estimate of the central black-hole 
mass can be gained from : Mj,^ ~ G -3 RBLnVg \ where Rblr 
is the radius of the BLR and V g is the velocity of the line- 
emitting gas, as traditionally estimated from the FWHM of 
the Hp emission line. 

Given the large number of physical processes which 
could potentially influence the dynamics of the BLR it is 
not immediately obvious that the gravitational potential of 
the central black-hole should be dominant (eg. Krolik 2001). 
However, several lines of evidence have recently shown that 
the dynamics of the BLR appear to be at least consistent 
with the virial assumption. Firstly, for the small number of 
objects for which it is possible to do so, Peterson & Wandel 
(2000) have shown that the motions of the broad-line gas are 
consistent with being virialized, with the velocity-widths of 
emission lines produced at different radii following the ex- 
pected V oc r -0 ' 5 relationship. Secondly, the black-hole mass 
estimates produced by the virial method are in good agree- 
ment with the predictions of the tight correlation between 
black-hole mass and stellar- velocity dispersion (Gebhardt et 
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al. 2000b; Ferrarese et al. 2001; Nelson et al. 2003; Green et 
al. 2003). Finally, using the virial estimator McLure & Dun- 
lop (2002) recently demonstrated, for 72 AGN at z < 0.5, 
that AGN host galaxies follow the same correlation between 
black-hole mass and bulge mass as local quiescent galaxies. 
Viewed in isolation, each of these lines of evidence might be 
seen as simply a consistency check. However, taken together 
they offer good supporting evidence for the validity of the 
virial assumption. 

Due to the fact that reverberation mapping measure- 
ments of Rblr are necessarily reliant on high-accuracy 
monitoring programmes lasting many years, at present such 
measurements are only available for a small sample of low- 
redshift AGN. Consequently, for estimating the black-hole 
masses for large samples of AGN, a crucial result arising 
from the Kaspi et al. (2000) study was that Rblr is strongly 
correlated with the AGN monochromatic continuum lumi- 
nosity at 5100 A . By exploiting this correlation it is there- 
fore possible to produce a virial black-hole mass estimate 
based purely on a luminosity and H/3 FWHM measurement. 
The availability of this technique has led to a proliferation 
of studies of AGN black-hole masses in the recent litera- 
ture. These studies have primarily focused on low-redshift 
(z < 0.5) AGN samples, investigating the relationships be- 
tween black-hole mass, the properties of the surrounding 
host galaxies and the spectral energy distribution of the 
central engine (e.g. Laor 1998, 2000, 2001; Wandel 1999; 
McLure & Dunlop 2001,2002; Lacy et al 2001; Dunlop et 
al 2003; Vestergaard 2004). However, the usefulness of the 
virial estimator based on the H/3 emission line is limited by 
the fact that H/3 is redshifted out of the optical at z > 0.8. 
Consequently, the use of this emission line to trace BLR ve- 
locities in high-redshift AGN requires infra-red spectroscopy, 
which is relatively observationally expensive and limited to 
the available atmospheric transmission windows. 

However, McLure & Jarvis (2002) and Vestergaard 
(2002) have recently demonstrated that this problem can be 
overcome by using the Mgll and Civ emission lines as rest- 
frame UV proxies for Hf3. In addition, McLure & Jarvis 
(2002) showed that Rblr is also strongly correlated with 
the monochromatic continuum emission at 3000 A, allow- 
ing black-hole mass estimates for high-redshift quasars to 
be made from a single optical spectrum covering the Mgll 
emission-line. 

The availability of the new rest-frame UV black-hole 
mass estimators has been exploited by recent studies to in- 
vestigate the black-hole masses of the most luminous quasars 
(Netzer 2003), the evolution of the black- hole mass - lumi- 
nosity relation (Corbett et al. 2003) and also to estimate the 
mass of the most distant known quasar at z = 6.41 (Willott, 
McLure & Jarvis 2003; Barth et al. 2003). In particular, 
using composite spectra generated from > 22000 2dF+6dF 
quasars, Corbett et al (2003) successfully demonstrated that 
the evolution of the black-hole mass - luminosity relation is 
too weak to explain the evolution of the z < 2.5 quasar 
luminosity function with a single population of long-lived 
objects. 

The recent publication of the SDSS first data release 
has provided publically available, fully calibrated, opti- 
cal spectra of > 17000 quasars in the redshift interval 
0.08 < z < 2.3. The flux calibration, long wavelength cov- 
erage (3800 A< A < 9200 A) and good spectral resolution 



(1.4 A/pix) of the SDSS spectra make them ideal for use in 
virial black-hole mass estimates. In this paper we study the 
black-hole mass distribution of quasars in the redshift range 
0.1 < z < 2.1 by utilizing continuum and emission-line fits 
to SDSS quasar spectra, together with calibrations of the 
virial mass estimators for H/3 and Mgll based on those de- 
rived by McLure & Jarvis (2002). Via cross-calibration with 
low-redshift AGN with actual reverberation mapping black- 
hole mass estimates (Kaspi et al. 2000) , it is known that the 
scatter in the virial mass estimators based on both the H(3 
and Mgll emission lines is a factor of 2.5 (see McLure & 
Jarvis 2002; appendix) and, given the potential uncertain- 
ties, it is probably unwise to trust individual virial black- hole 
mass estimates to within a factor of ~ 3. However, because 
there is no systematic off-set between the virial estimators 
and the reverberation mapping results, the virial mass esti- 
mators are expected to produce statistically accurate results 
when applied to large samples such as the SDSS. 

The structure of the paper is as follows. In Section 2 we 
described the criteria used for selecting our sample from the 
SDSS Quasar Catalogue II of Schneider et al. (2003). In Sec- 
tion 3 the basic distribution of the SDSS quasar black-hole 
masses as a function of redshift is presented. This new infor- 
mation is then used to address important questions regard- 
ing the virial black-hole mass estimator and the the accretion 
rate of luminous quasars. Firstly, a recent study by Netzer 
(2003) based on the Civ and H/3 virial mass estimators has 
suggested that many luminous quasars at z > 2 may har- 
bour central black holes with masses in excess of 10 10 M . 
As highlighted by Netzer (2003), since galaxies with suit- 
ably large velocity dispersions (<r>600 km s _1 ) are not ob- 
served at low redshift, this result raises questions about ei- 
ther the reliability of the virial mass estimator, or the form 
of the bulge:black-hole mass relation, at high redshift. We 
re-address this question using our new SDSS quasar black- 
hole mass estimates based on the Mgll calibration of the 
virial mass estimator. 

Secondly, recent studies have shown that at low-redshift 
(z < 0.5), luminous broad-line quasars are accreting be- 
low the Eddington limit (Floyd et al. 2003; Dunlop et al. 
2003; McLeod & McLeod 2001). However, it has been sug- 
gested (eg. Mathur 2000) that narrow-line Seyfert galaxies 
(NLS1) are the low-redshift analogs of the early evolution of 
luminous quasars, where the central black hole is still in an 
exponential growth phase, accreting beyond the Eddington 
limit. Consequently, if NLS1 are the low-redshift analogs of 
the early evolution of luminous quasars, then at high redshift 
we might expect to see luminous quasars accreting above the 
Eddington limit. Using the black-hole mass and bolometric 
luminosity estimates for our SDSS sample, we investigate 
whether the Eddington luminosity is still a physically rele- 
vant limit to quasar accretion rates at z < 2. 

In Section 4 we re-investigate the form of the local 
black-hole mass function within the context of the latest 
results on the Mbh — Li, u i ge and Mt,h — a correlations. 
In Section 5 this information is combined with the new 
SDSS quasar black-hole masses and the 2dF quasar lumi- 
nosity function to estimate the activation fraction of z ~ 2 
black holes as a function of mass, and deduce lower limits 
on the lifetimes of the most luminous quasars. In Section 
6 we summarize our main results and conclusions. Full de- 
tails of the calibration of the virial black-hole mass estima- 



© 0000 RAS, MNRAS 000, 000-000 



The cosmological evolution of quasar black-hole masses 3 




z 

Figure 1. Virial black-hole mass estimate versus redshift for our full SDSS quasar sample. Broad-line quasars (FWHM > 2000 km s _1 ) 
are shown as grey symbols, while narrow-line objects (FWHM< 2000 km s _1 ) arc shown as black symbols. The mean black-hole masses 
within Az = 0.1 bins arc shown as filled circles (standard errors are smaller than the symbols except for lowest redshift bin). The vertical 
dotted line highlights the switch from using the H/3-based to the Mgll-based virial mass estimator at z = 0.7. The horizontal solid line 
marks a black-hole mass of 3 X 10 9 Mq, the maximum mass observed at low redshift (see text for discussion). 



tors and the emission-line modelling of the SDSS quasar 
spectra are provided in the Appendix. Unless otherwise 
stated all calculations assume the following cosmology : 
H = 70 km s^Mpc" 1 , Q m = 0.3, Q A = 0.7. 



2 THE SAMPLE 

The sample of quasars analysed in this paper is drawn from 
the SDSS Quasar Catalog II (Schneider et al. 2003). This 
catalog is based on the SDSS first data release, and features 
16713 quasars with Mi(AB) < —22 in the redshift interval 
0.08 < z < 5.41. The first stage in constructing our sam- 
ple was to select all objects from the Schneider et al. catalog 
in the redshift range 0.08 < z < 2.1, with the upper red- 
shift limit imposed to ensure sufficient wavelength coverage 
to reliably determine the Mgll emission-line widths. The re- 
sulting list of 14181 quasars was then passed through the 
emission-line modelling software described in the Appendix. 

The final sample used in the analysis consists of 12698 
quasars, or 90% of the z < 2.1 quasars in the Schneider et 
al. catalogue. The remaining 10% of objects were excluded 
for having unreliable modelling results due to low signal-to- 
noise, particularly weak emission lines or having been af- 
fected by artifacts in the spectra. In conclusion, although 
the SDSS Quasar Catalog II is not statistically complete 
(Schneider et al. 2003), our final sample of 12698 quasars is 
clearly representative of optically luminous quasars in the 
redshift interval 0.1 < z < 2.1. 



3 THE EVOLUTION OF QUASAR 
BLACK-HOLE MASSES 

In Fig 1 we show the distribution of virial black-hole mass 
estimates versus redshift for our full SDSS quasar sample. 
Broad-line quasars (either H/3 or Mgll emission line with 
FWHM> 2000 km s _1 ) are shown as grey symbols and 
narrow- line objects (FWHM< 2000 km s _1 ) are shown as 
black symbols. The mean black-hole masses within Az — 0.1 
bins are shown as filled circles, and the vertical dotted line 
at z — 0.7 marks the transition between the use of the H[3 
and Mgll based virial estimators. There can be seen to be a 
smooth transition between the two virial mass calibrations, 
and more details concerning the quality of the agreement 
between the H/3 and Mgll based virial estimators are pro- 
vided in the Appendix. Finally, the solid horizontal line is at 
a black-hole mass of 3 x 10 9 M Q , which is representative of 
the maximum mass which has been dynamically measured 
within the local galaxy population (Ford et al. 1994; Tad- 
hunter et al. 2003). 



3.1 Black- hole mass as a function of redshift 

The mean black-hole mass can be seen to increase as a 
function of redshift, rising from ~ 10 s M Q at z ~ 0.2 to 
~ 10 9 Mq at z ~ 2.0. However, this is entirely as ex- 
pected due to a combination of the effective flux-limit of 
the sample and the role of the quasar continuum luminosity 
within the virial mass estimator. In the event that the mean 
FWHM remains approximately constant with redshift, the 
mean virial black-hole mass estimate is required to increase 
as < \L\ > ' 6 by construction (Eqn A6 & Eqn A7). The 
increase in mean black-hole mass with redshift seen in Fig 
1 is entirely consistent with this effect. 
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Figure 2. Panel A shows bolometric luminosity as a fraction of the Eddington limit versus redshift for the full sample. The grey symbols 
are broad-line quasars (FWHM > 2000 km s -1 ) and the black symbols arc narrow-line objects (FWHM< 2000 km s^ 1 ). The mean 
values of L bo i/L Edd within Az = 0.1 bins are shown as filled circles. The solid horizontal line marks the Eddington limit. Panel B is 
identical to panel A except that only the broad-line quasars have been plotted. 



3.2 Evidence for a limiting black-hole mass 

Amongst the sample of local galaxies for which the central 
black-hole mass has been measured via gas or stellar dynam- 
ics, the most massive black holes which have been found to 
date are in M87 (Ford et al. 1994) and Cyg A (Tadhunter 
et al. 2003), both with M bh ~ 3 x 10 9 M . 

Given the relatively small cosmological volume from 
which this sample of local galaxies is drawn, it is of interest 
to investigate whether ~ 3 x 10 9 M represents a genuine 
limiting mass. Furthermore, the Mbh — a and Mbh — Lbuige 
correlations, defined by the local galaxies with dynamical 
black-hole mass measurements, make definite predictions for 
what the maximum black-hole mass is likely to be. For exam- 
ple, the distribution of stellar-velocity dispersions for early- 
type galaxies from the SDSS (Bernardi et al. 2003; Seth et 
al. 2003) show that velocity dispersions of >400 km s _1 are 
exceptionally rare among the local galaxy population. Com- 
bined with the latest determination of the M b h — cs relation 
(Tremaine et al. 2002), this implies a limiting black- hole 
mass of ~ 2 x 10 9 M© , in excellent agreement with the largest 
black-hole masses measured to date. Consequently, if the 
virial black-hole mass estimator predicts that a large pop- 
ulation of luminous quasars at high redshift harbour black 
holes with masses > 10 10 Mq , the implied velocity disper- 
sions of a ~ 600 km s~ would be in direct conflict with the 
known properties of early-type galaxies. 

This issue was recently addressed by Netzer (2003), who 
compiled virial black-hole mass estimates (based on H/3 and 
Civ emission-line widths) for a heterogeneous sample of 724 
quasars in the redshift interval < z < 3.5. The results 
of this study indicated that a significant number of lumi- 
nous (ALi35o > 10 46 ' 5 ergs s _1 ) quasars are powered by black 
holes with masses >10 Mq. Consequently, Netzer (2003) 
concluded that either the virial mass estimator (principally 
the Rblr — ALa correlation) was invalid at z>2 or, alter- 
natively, that the relationship between black-hole and bulge 
mass has a different form at high redshift, presumably be- 
cause the quasar hosts are still in the process of forming. 

However, the new results presented here for the SDSS 
quasars do not appear to support this conclusion. It can be 



seen from Fig 1 that the upper envelope of black-hole masses 
for the SDSS quasars is fully consistent with a maximum of 
~ 3 x 10 9 M© , with a negligible number of quasars exceeding 
this limit at z < 1.5. At the high-redshift end of the sample 
(1.5 < z < 2.1) there are only a small number of quasars 
with virial black-hole mass estimates > 3 x 1O 9 M0, less 
than 5% of the SDSS quasars with z > 1.5, and only two 
objects are found to have virial black-hole mass estimates 
in excess of 10 10 Mq. When the level of scatter associated 
with the virial mass estimator (~ 0.4 dex) is taken into ac- 
count, it is clear that there is no conflict between the SDSS 
quasar black-hole masses and either the locally observed 
maximum black-hole mass or, perhaps more significantly, 
the maximum black-hole mass predicted from combining the 
stellar-velocity dispersion distribution of early-type galaxies 
(Bernardi et al. 2003) with the Mbh — a relation. 

The reason for the discrepancy between the new results 
presented here and those of Netzer (2003), hereafter N03, 
is not entirely clear. The distribution of black-hole mass as 
a function of redshift determined by N03 shows an upper 
limit consistent with Mbh — 10 M© , with a small number 
of objects with estimated masses exceeding this limit. In 
contrast, the results for the SDSS quasars presented in Fig 
1 are consistent with a limit of ~ 3 x 10 9 Mq , a difference of 
0.5 dex. The discrepancy cannot be due to the use of differ- 
ent calibrations of the Rblr — \L\ correlation. If fact, if our 
virial mass estimator is adjusted to the same Rblr — \L\ 
slope (AL 9 ) and the same geometric normalization (N03 
adopt / = v3/2, we adopt f = 1, see Eqn Al) our black- 
hole mass estimates for the most luminous SDSS quasars 
are actually reduced by ~ 0.2 dex. Presumably the cause 
of the discrepancy must therefore lie either in the different 
methods used to determine the quasar continuum luminosi- 
ties, or N03's use of Civ as a proxy for H/3 instead of the 
Mgll FWHM measurements adopted here. In this context 
we note that Vestergaard (2004) recently found evidence for 
a similar limiting black-hole mass to that found here, using a 
smaller sample of 2 < z < 4 quasars and the Civ-based virial 
mass estimator. 

In conclusion, based on the H(3 and Mgll virial black- 
hole mass estimates for > 12000 SDSS quasars presented 
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Figure 3. The distribution of the full SDSS sample on the accre- 
tion rate — black-hole mass plane. The grey symbols are broad- 
line quasars (FWHM > 2000 km s _1 ) and the black symbols are 
narrow- line objects (FWHM< 2000 km s" 1 ). The solid line shows 
the location of the Eddington limit (Lboi = L^dd)) while the 
dashed line highlights an accretion rate of ten solar masses per 
year, below which 94% of the sample lie. 

here, we find no evidence that z < 2 quasars harbour cen- 
tral black-holes with masses in excess of the maximum of 
~ 3 x 10 Mq observed in the local Universe. This consis- 
tency between the most massive black-holes at z ~ 2 and the 
most massive black-holes at z ~ is important given that 
it is known from low-redshift studies that the host galax- 
ies of these z ~ 2 quasars must be, or at least evolve into, 
massive early-type galaxies (Dunlop et al. 2003; McLure & 
Dunlop 2002; McLeod & McLeod 2001). Reassuringly, the 
black-hole mass distribution of the SDSS quasars presented 
here is entirely consistent with their host galaxies having 
similar properties to those of low-redshift early-types; i.e. 
a<400 km s" 1 and L<10L*. 



3.3 Quasar accretion rates 

In panel A of Fig 2 we plot the distribution of quasar Edding- 
ton ratios (L bo i/L E dd) as a function of redshift for the full 
sample (narrow + broad-line objects). The mean Lboi/LEdd 
values within Az — 0.1 bins are shown as the filled cir- 
cles, with the solid horizontal line highlighting the Edding- 
ton limit (Lboi/LEdd = 1). Panel B is identical to panel 
A except that only the broad-line objects are plotted. Sev- 
eral features of the distribution of the SDSS quasars on the 
Lboi/LEdd — z plane are worthy of individual comment and 
are discussed in turn below. 

In common with the distribution of black-hole masses, 
it can be seen that the mean Eddington ratio of the SDSS 
quasars increases with redshift, rising from Lboi — 0.15LEdd 
at z ~ 0.2 to L b oi — 0.5L E dd at z ~ 2.0 (details of how 
bolometric luminosities have been estimated from the quasar 
continuum luminosities are provided in the Appendix) . How- 
ever, due to the fact that the Eddington luminosity is di- 
rectly proportional to the estimated black-hole mass, the 
mean Eddington ratio is required to scale as Lboi/LEdd oc 



Lboi A in the event that the mean FWHM is approximately 
constant with redshift. The increase in the mean Eddington 
ratio with redshift seen in Fig 2 is entirely consistent with 
this expectation. 

The distribution of Eddington ratios displayed by the 
SDSS quasars clearly shows that assuming that all lumi- 
nous quasars are accreting at their Eddington limit is a poor 
approximation. This result is important because it is often 
assumed that optically luminous quasars are accreting at 
their Eddington limit within models of quasar evolution. 
In fact, our results show that a constant accretion rate of 
~ 0.3 — > 0AL E dd is a much better approximation than using 
the Eddington luminosity. 

In panel A of Fig 2 it can be seen that the majority 
of the narrow-line objects (black symbols) are estimated 
to be accreting at super-Eddington rates. This result is in 
agreement with numerous studies of low-redshift narrow- 
line Seyfert 1 galaxies (NLS1) which, in the most widely 
accepted model, are thought to be examples of low black- 
hole mass/high accretion rate AGN. However, it is also im- 
portant to remember that the narrow-line portion of the 
SDSS quasar sample is likely to be an admixture of genuine 
narrow-line objects (i.e. low black-hole mass/high accretion 
rate) and intrinsically broad-line quasars which appear as 
narrow-line objects simply because they are viewed close to 
pole-on. Indeed, as discussed in the appendix, if there is 
any significant disk-component to the FWHM of the broad 
emission lines then some level of contamination by pole-on 
sources is inevitable. Within this context it is noteworthy 
that we do not see any trend for an increase in the num- 
bers of narrow-line objects with redshift, as might be ex- 
pected under the hypothesis that NLSls are the low-redshift 
analogs of the early evolutionary stages of luminous broad- 
line quasars (eg. Mathur 2000). On the contrary, the fraction 
of narrow-line objects in our SDSS sample decreases from 
~ 10% at z ~ 0.2 to only ~ 1% at z ~ 2 (see Fig 1). In fact, 
this trend is at least qualitatively consistent with the pre- 
dictions of a simplified model of the BLR which depends on 
both orientation and luminosity. Firstly, if it is assumed that 
all quasars are intrinsically broad-line objects, and that the 
BLR has a flattened geometry, then narrow-line objects are 
simply those which are viewed close to pole-on (< 15°). In a 
flux-limited sample the proportional of narrow-line objects 
is then predicted to fall with redshift if the opening angle of 
the obscuring torus is dependent on the AGN luminosity; i.e 
8 oc L 0,5 (the so-called "receding torus" model e.g. Simpson 
1998). In this simplified scheme the redshift evolution of the 
narrow-line fraction seen in Fig 1 is roughly consistent with 
a receding torus model with a half-opening angle of ~ 45° at 
the low-redshift/luminosity end of the SDSS sample, rising 
to ~ 80° at the high-redshift/luminosity end. 

It can be seen from panel B of Fig 2 that at redshifts 
1 < z < 2 the upper envelope of the Eddington ratios for 
the broad- line quasars is virtually flat at Lboi/LEdd — 1- 
Indeed, the fact that we do not see significant numbers of 
quasars which are accreting at greater than the Eddington 
limit, with virtually none accreting at > 2L E dd, can be taken 
as strong evidence that the Eddington limit is a physically 
meaningful upper limit to the accretion rate in broad-line 
quasars, at least at z < 2 . 

In Fig 3 we show the distribution of the SDSS quasars 
on the accretion rate — black-hole mass plane, where the 
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Figure 4. Two independent estimates of the redshift zero dor- 
mant black-hole mass function for early-type galaxies. The solid 
line is a transformation of the if— band early-type galaxy lumi- 
nosity function of Kochanck ct al. (2001) using the Mbh — Lbulge 
correlation of McLure & Dunlop (2002). The dashed line is a 
transformation of the SDSS early-type galaxy stellar-velocity dis- 
persion function of Seth et al. (2003) using the Trcmainc ct al. 
(2002) fit to the M bh - a correlation. 

accretion rates have been calculated assuming a canonical 
mass-to-energy conversion efficiency of e = 0.1. The solid 
line in Fig 3 shows the location of the Eddington limit 
(Lboi = L E dd), while the dashed line highlights an accretion 
rate of ten solar masses per year. It can immediately be seen 
from Fig 3 that the vast majority (94%) of the SDSS quasars 
are accreting material at rates of M < M < 10 M per 
year. However, perhaps the most striking feature of Fig 3 
is the near total lack of quasars with black-hole masses of 
Mbh > 10 9 Mq and accretion rates close to the Eddington 
limit. This feature must be regarded as significant because 
it is extremely unlikely that large numbers of such objects 
are missing from the SDSS quasar sample. We note here that 
this situation is at least qualitatively consistent with models 
of black-hole growth in which the exponential, Eddington- 
limited, growth of the black-hole is terminated by a physical 
limit to the amount of material which can be supplied for ac- 
cretion (eg. Archibald et al. 2002; Granato et al. 2003, 2001). 
In models of this type, following the end of the exponential- 
growth phase the black hole is free to remain active as a 
luminous quasar for ~ 10 8 years, accreting below the Ed- 
dington limit, without leading to the production of large 
numbers of black holes with masses > 10 10 Mq which are 
not observed locally. This type of scenario is also consistent 
with the results shown in Fig 1 & Fig 2, and is discussed in 
more detail in Section 5. 



4 THE LOCAL MASS FUNCTION OF 
DORMANT BLACK HOLES 

Following the discovery that supermassive black holes ap- 
pear to be ubiquitous at the centres of massive local galax- 
ies, it has been of interest to use this information to estimate 



the mass function of dormant black holes in the local Uni- 
verse. From a theoretical perspective it was shown by Soltan 
(1982) that the local mass density of dormant black holes 
could be calculated from the total radiated energy of opti- 
cal quasars, which in turn can be calculated from the quasar 
source-count distribution. Traditionally such estimates pre- 
dict local black-hole mass densities a factor of a few lower 
(Soltan 1982; Chokshi & Turner 1992; Yu & Tremaine 2002) 
than those estimated from the X-ray background (eg. Fabian 
& Iwasawa 1999), suggesting that the dominant fraction of 
the quasar population is optically obscured. Given that a 
principal objective of quasar evolution models is to simulta- 
neously explain the evolution of both the quasar luminosity 
and black-hole mass functions, it is clearly important to have 
a robust measurement of the form of the local mass function 
of dormant black holes. 

Previous estimates of the local black-hole mass function 
in the literature have relied on either the Mbh — Lbulge corre- 
lation to make the transformation from the local galaxy lu- 
minosity function (e.g. Salucci et al. 1999), or the correlation 
between radio luminosity and black-hole mass observed in 
local early- type galaxies (e.g. Franceschini et al. 1998). More 
recently, both Yu & Tremaine (2002) and Aller & Richstonc 
(2002) exploited the tight Mbh — & correlation to derive an 
estimate of the local mass density of dormant black holes. 
While the Yu & Tremaine (2002) calculation was based di- 
rectly on the SDSS stellar-velocity dispersion function, the 
Aller & Richstone (2002) calculation used the Faber-Jackson 
relation (Faber & Jackson 1976) to first estimate the dis- 
persion function from the galaxy luminosity function. De- 
spite this difference in approach, both studies found consis- 
tent results for the local black-hole mass density (early+late 
types) with estimates of (2.4 ± 0.8) x 10 5 M Q Mpc~ 3 and 
(2.9 ± 0.5) x 10 5 MgMpc" 3 for Aller & Richstone and Yu & 
Tremaine respectively (Ho = 70 km s _1 Mpc _1 ). 

In Section 5 we will proceed to estimate the activation 
fraction of supermassive black holes at z ~ 2. This calcula- 
tion requires a knowledge of the actual functional form of the 
local dormant black-hole mass function, information which 
is unavailable from the previous studies in the literature. At 
z ~ 2 we are exclusively interested in the form of high-mass 
end of the local black-hole mass function, and consequently, 
in this section we use the Mbh — a and Mbh — L bu ige relations 
to derive two independent estimates of the local black-hole 
mass function for early-type galaxies. 

4.1 The velocity dispersion estimate 

Seth et al. (2003) recently performed a detailed analysis of 
the stellar-velocity dispersion function of 9000 early-type 
galaxies drawn from the SDSS. As part of this analysis Seth 
et al. provide a fitting formula for the dispersion function, 
similar to a Schechter function in form, which is governed 
by four free parameters (</>*, ct*, a, f3). As the starting point 
for our calculation we adopt the best fit determined by Seth 
et al. to the observed dispersion function, including mea- 
surement errors, which is described by the parameter values 
(0.002 Mpc~ 3 , 88 km s"\ 6.5, 1.8). Secondly, we then apply 
a change of variables to convert the velocity dispersion func- 
tion into a black-hole mass function using the latest version 
of the M bh - a correlation (M bh oc 0- 4 O2 ± - 32 ) as derived 
by Tremaine et al. (2002). The final step in the process is 
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to convolve this estimate of the black-hole mass function 
with a gaussian of a — 0.27 dex to account for the associ- 
ated scatter of the Mbh — a relation. The estimate of the 
local black-hole mass function produced via this method is 
shown as the dashed line in Fig 4. 

4.2 The bulge luminosity estimate 

Our second approach involves a conversion of the local 
galaxy luminosity function into an estimate of the local 
black-hole mass function using the Mbh — Lbuige correla- 
tion. This approach was previously used by Salucci et al. 
(1999), who derived an estimate of the local black-hole mass 
function via a two-step process. The first stage was the con- 
struction of the local bulge/spheroid mass function from rep- 
resentative optical galaxy luminosity functions, segregated 
by morphological type, using an appropriate mass-to-light 
ratio and an estimate for each morphological type of the 
bulge:total light ratio. The second stage was to then convert 
the spheroid mass function into the black-hole mass function 
by convolution with a log gaussian Mbh /Mbuige distribution, 
centred on < Mbh /Mbuige >= —2.6 with a — 0.3 dex. The 
Salucci et al. (1999) study demonstrated that at the high- 
mass end (Mbh > 10 8 M Q ), the local dormant black-hole 
mass function should be entirely dominated by the contri- 
bution from early-type (E+SO) galaxies. 

We are now in a position to re-visit this calculation 
using improved information on the luminosity function of 
early-type galaxies and the on form of the black-hole mass 
- bulge mass relation. For the purposes of our estimate we 
adopt the 2MASS early-type if— band luminosity function 
of Kochanek et al. (2001). The Kochanek et al. early-type 
luminosity function is conveniently segregated to included 
only SO and elliptical galaxies, making it well matched to 
the SDSS velocity dispersion function of Seth et al. (2003). 
Moreover, due to the if— band selection, we can be confi- 
dent that the galaxy luminosities trace the mass of the old 
stellar spheroid population, which has been shown recently 
to correlate with central black-hole mass with a comparable 
scatter to the Mbh — o correlation; i.e ~ 0.3 dex (McLure & 
Dunlop 2002; Erwin et al. 2003; Marconi & Hunt 2003). 

Specifically, we adopt the K~ band luminosity function 
parameters derived by Kochanek et al. (2001) for their early- 
type (E+S0) sample: a = -0.92, <j> = 0.45 x 10" 2 /i? 00 Mpc~ 3 
and Mk = —23.53 + 51ogftioo- After conversion of the lu- 
minosity function parameters to our adopted cosmology, the 
value of M* K was further correct by Am = —0.2 to trans- 
form from isophotal to total magnitudes (Kochanek et al. 
2001). In order to make the transformation from the lumi- 
nosity function to the black-hole mass function we adopt 
the best-fit to the Mbh — Lbuige correlation from McLure 
& Dunlop (2002), which was determined from a sample of 
exclusively elliptical galaxies with reliable dynamical black- 
hole mass measurements. After accounting for the change in 
cosmology, this relation becomes: 

log M bh = 1.25(±0.05) log ^ - 5.76(±0.53) (1) 

where an average colour of R — K — 2.7 for early-type E/S0 
galaxies has been adopted (GISSEL98 models; Bruzual & 
Chariot 1993), and we have also assumed for this simplified 
calculation that the K— band bulge:total luminosity ratio for 



both elliptical and SO galaxies is unity. Following convolution 
with a gaussian with a — 0.3 dex to account for the scatter 
around the M b h — L bu i ge correlation, the final estimate of 
the early-type dormant black-hole mass function is shown 
as the solid line in Fig 4. 



4.3 Integrated black-hole mass density 

It can be seen from Fig 4 that the two independent estimates 
of the local early-type dormant black-hole mass function are 
in good agreement at the high- mass end (Mbh > 10 s Mq). 
The turn-over of the velocity-dispersion based estimate at 
low masses is an artifact of the SDSS velocity dispersion 
function, which turns over at a < 125 km s _1 due to the faint 
magnitude limit of the SDSS early-type sample (Bernardi et 
al. 2003). Finally, if we integrate our estimate of the local 
black-hole mass function based on the K — band luminosity 
function, our estimate of the local black-hole mass density 
in early type galaxies is: 

Pbh = (2.8 ± 0.4) x 10 5 MgMpc -3 (2) 

which, considering the potential uncertainties and simplified 
nature of this calculation, can be seen to be in good agree- 
ment with the values determined by both Yu & Tremaine 
(2002) and Aller & Richstone (2002). 

The good agreement between the two estimates of the 
local dormant black-hole mass function found here is appar- 
ently in contrast to one of the findings of Yu & Tremaine 
(2002), who found that using the Mbh — Lbuige relation led to 
a black-hole mass density estimate more than a factor of two 
greater than their velocity dispersion estimate. Although 
small differences are introduced by our different choice of 
Mbh — Lbuige relation and galaxy luminosity function, wc 
note here that the principal source of this apparent discrep- 
ancy is that we have adopted virtually identical levels of 
scatter in the M b h — a and Mbh — Lbuige relations (~ 0.3 dex). 
In contrast, Yu & Tremaine adopted a scatter of 0.5 dex in 
the Mbh — Lbuige relation, based on the fit to B— band data 
by Kormendy & Gebhardt (2001). If the Yu & Tremaine 
Mbh — Lbuige density estimate is adjusted to reflect a scatter 
of ~ 0.3 dex then their two density estimates become consis- 
tent, differing by less than 50%. Although the relative levels 
of scatter associated with the Mbh — a and Mbh — Lbuige 
relations is currently a question of some debate, our choice 
of adopting a low scatter of ~ 0.3 dex in the Mbh — Lbuige 
relation is justified, especially for early-types, by the recent 
results of several detailed studies (McLure & Dunlop 2002; 
Erwin et al. 2003; Marconi & Hunt 2003). 



5 CONSTRAINING THE ACTIVATION 

FRACTION OF SUPERMASSIVE BLACK 
HOLES 

Although the completeness of the SDSS Quasar Catalog II 
has not been fully quantified (Schneider et al. 2003), there is 
no obvious reason to suspect that the SDSS quasars are not 
at least representative of luminous optical quasars. There- 
fore, it is possible to use our new SDSS quasar black-hole 
mass estimates to study the distribution of quasar black-hole 
masses as a function of luminosity. Furthermore, by using 
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Figure 5. Each panel shows the distribution of absolute B— band magnitude for the 1069 quasars from our SDSS sample in the redshift 
interval 1.90 < z < 2.1. The black histograms show the distribution of sub-sets of these quasars whose black- holes are more massive 
than the thresholds shown in the top-left corner of each panel. 



the optical quasar luminosity function to calculate the num- 
ber density of quasars brighter than a given absolute mag- 
nitude, it is then possible to estimate the number density of 
active black-holes as a function of mass. In this section we 
perform this calculation using our black-hole mass estimates 
for the SDSS quasars at z ~ 2. 

For the purposes of this calculation we adopt the quasar 
luminosity function as determined by Boyle et al. (2000) 
for their combined 2dF+LBQS sample of more than 6000 
quasars in the redshift range 0.35 < z < 2.30. The Boyle 
et al. luminosity function is defined in terms of absolute 
B— band magnitudes, so we initially convert the bolomet- 
ric luminosity estimates for our SDSS quasars into abso- 
lute B— band magnitudes, using a transformation defined 
by > 300 quasars common to both the SDSS and the 2dF 
10K quasar catalog (Croom et al. 2001; see Appendix). 

In Fig 5 we show the distribution of absolute B— band 
magnitude for 1069 SDSS quasars in the redshift interval 
1.90 < z < 2.10. The filled histograms in each panel 
show the distribution of absolute B— band magnitudes for 
three sub-sets of these objects which have black-hole masses 
greater than our chosen mass thresholds of 10 8 ' 5 M Q , 10 9 M 
& 10 9 ' 5 M Q . The solid line at M B = -26 shows the location 
of the mi(AB) = 19.1 (Schneider et al. 2003) flux limit of 
the main SDSS multi-colour quasar selection algorithm at 
z — 2.1. Provided that at a given quasar luminosity the 
SDSS is not significantly biased towards including only the 
most massive black-holes, the ratio of the filled to un-filled 
histograms in Fig 5 provides a method for determining lower 
limits on the number densities of active black holes at z ~ 2. 

We perform this calculation in two steps. Firstly, for 
each of the three mass thresholds we calculate the fraction 
/ of quasars with black-hole masses above the threshold, in 
three magnitude ranges: —25 < Mb < —26, —26 < Mb < 
-27 & -27 < M B < -28.3 (switching to 0.5 magnitude 
bins makes little difference to the final estimate). The lower 
limit of Mb = — 25 has been chosen to conservatively ensure 
that we still have sufficient numbers of quasars in the filled 
histograms to determine / in each magnitude bin to better 
than 10% accuracy. Although Mb = — 25 is one magnitude 
fainter than the SDSS quasar selection algorithm flux limit 



at z ~ 2, we are making the explicit assumption that the 
ratio of the filled to un-filled histograms would remain un- 
changed, even in the event that the SDSS quasar selection 
was complete to Mb = —25. The second step in the process 
is to use the Boyle et al. (2000) quasar luminosity function 
to calculate the number density of quasars {pact) in the three 
magnitude bins at z = 2. For this calculation we adopt the 
Boyle et al. polynomial fit to the luminosity function, con- 
verting the absolute magnitude limits and number densities 
to our adopted cosmology. Finally, the estimated number 
densities of active black-holes for each of the mass thresh- 
olds is simply the sum of / x p ac t over the three magnitude 
bins. 

The estimates of the number densities of active super- 
massive black holes at z ~ 2 using the above method are 
shown in panel A of Fig 6, for the three black-hole mass 
thresholds. Also shown are the z ~ cumulative black-hole 
mass functions based on our M^h — L ou i ge and M„h — a es- 
timates from Section 4. This plot indicates that the fraction 
of present-day black holes which are in place, active, and op- 
tically unobscured at z ~ 2 is ~ 0.005 for M bh > 1O 8 ' 5 M , 
~ 0.02 for M bh > 1O 9 M , and ~ 0.05 for M bh > 1O 9 ' 5 M . 

How should we interpret these ratios? The actual ob- 
served number density ratio at a given mass threshold, R, 
must be a product of three physical ratios, namely: 



R = {n z =2/n z= o) x (tq/t epoch ) x (n vis /n obs ) 



(3) 



where n z =2 is the number density of black holes (above the 
appropriate mass threshold) already in place at z = 2, n z =o 
is the corresponding number density at the present day, ig is 
the net lifetime of a black hole (above the appropriate mass 
threshold) in the optically luminous quasar phase, t epoc h is 
the total duration of the luminous quasar epoch, and the 
ratio {n v is / n oba ) describes the effect of geometrical obscura- 
tion by dusty tori (note the effect of a completely obscured 
phase is different, and is implicitly taken care of by the def- 
inition of to). 

Thus, while at first sight it might appear that panel 
A in Fig 6 indicates that the fraction of present-day black 
holes already in place at z ~ 2 is an increasing function 
of mass, this could equally well reflect a mass-dependence 
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Figure 6. In each panel the solid and dashed lines show the cumulative local dormant black-hole mass functions as derived in Section 4 
using the — L bu i ge and — a relations respectively. In panel A we show the estimated number densities of active supermassivc 
black-holes at z ~ 2 for three mass thresholds (> 10 s Mq, > 10 9 Mq & > 10 9 ' 5 M ). In panel B these number densities have be adjusted 
to account for a possible relationship between mean quasar lifetime and black-hole mass (see text for discussion). In Panel C the number 
densities of supermassivc black-holes have been boosted by a factor of two to conservatively account for geometric obscuration. 



of either net quasar lifetime, or torus opening angle. Some 
interesting limiting values can, however, be deduced with 
reasonable confidence. For example, adopting t epoc h = 2 Gyr 
(corresponding to the redshift range 1.5 < z < 3.5, the full 
width half maximum of the peak in luminous quasar number 
density - see, e.g. Richstone et al. (1998); Osmer 2003), we 
can infer a lower limit on the optical quasar lifetime of the 
most massive black holes M > 1O 9 ' 5 M of t Q > 1 x 10 8 
yr if there are no type II quasars, or tQ > 2 x 10 8 yr if 
(n vis /n obs ) ~ 0.5 as found at the bright end of the radio- 
loud quasar population. Such numbers seem very reasonable, 
agreeing well with, for example, the recent findings of Yu 
& Tremaine (2002). The fact that t Q > 1 x 10 8 yr differs 
spectacularly from, for example, the conclusion reached by 
Richstone et al. (1998) that tQ ~ 1 x 10 6 yr, can be traced 
to the fact that our inferred quasar black-hole masses are on 
average larger (because we have not been forced to assume 
Eddington limited accretion) and the fact our analysis is 
confined to the very high-mass end of the black-hole mass 
function. 

In fact, as a result of their analysis, Yu & Tremaine 
(2002) find evidence for a mass dependence of tQ which can 
largely explain the apparent mass dependence of activity ra- 
tio shown in panel A of Fig 6. Specifically, Yu & Tremaine 
determined that the mean quasar lifetime should be an in- 
creasing function of black-hole mass, rising from 3 x 10 7 yr 
at 10 8 M o to 5 x 10 8 yr at 1O 9 ' 5 M . For illustrative pur- 
poses, we show in panel B of Fig 6 the effect of correcting 
the estimated number densities at z ~ 2, using the appro- 
priate values of tQ as a function of mass (from Fig. 5 of Yu 
& Tremaine) and again adopting t epoc h = 2 Gyr. This has 
the effect of removing all significant mass dependence, and 
raises the apparent black-hole activation ratio to > 0.2 for 
all three mass thresholds. Applying a conservative correction 
of ~ 2 for geometric obscuration raises this ratio to > 0.4 
(see panel C of Fig 6). Given that our z ~ 2 census is un- 
doubtedly incomplete (being confined to the very luminous 
end of the quasar population) this lower limit is sufficiently 
close to unity to imply that the majority of supermassive 



black holes are in place by z ~ 2, and to leave little, if any, 
room for an additional substantial population of completely 
obscured growing black holes at the high-mass end. This 
finding is in accord with the conclusions of Yu & Tremaine 
(2002), and Fabian (2003) that the apparently substantial 
obscured population of AGN required to explain the X-ray 
background is to be found at lower black hole masses than 
those probed by the present study. 

Finally, we note that values of quasar lifetime (£q) as 
large as a significant fraction of a Gyr as discussed above, 
have often been rejected as unfeasible for the most mas- 
sive black holes M > 10 9 Mq . This is because, if accretion 
is assumed to proceed at, or close to, the Eddington limit, 
extreme black hole masses M > 10 10 Mq will then be pro- 
duced in only 1 x 10 8 yr, and as discussed in Section 3.2 
we find no evidence for any significant population of such 
extreme-mass black holes within the SDSS quasar sample. 
However, the lack of such extreme objects can be reconciled 
with our inferred lower limit on quasar lifetime tQ > 2 x 10 8 
yr by revisiting the results shown in Fig 3. Here it can be 
seen that continued Eddington-limited exponential growth 
towards M ~ 10 10 M requires mass accretion rates which 
rapidly approach lOOMoyr -1 . In contrast, allowing for the 
known scatter in the calculation of bolometric luminosities 
(Elvis et al. 1994), Fig 3 provides good evidence that the 
most massive black holes are not accreting matter at a rate 
significantly in excess of lOMeyr" 1 . Clearly, a 10 9 M black 
hole can continue to accrete matter at this rate, and pro- 
duce bright quasar light for a substantial fraction of a Gyr 
without producing a final black-hole mass in excess of a few 
times 1O 9 M . 

In conclusion, our results imply that the majority of 
the most massive black-holes are in place by z ~ 2 and that 
the most massive black holes accrete the bulk of their final 
mass (i.e. as they grow from ~ 10 8 ' 5 M to ~ 10 9 ' 5 M ) as 
optically luminous quasars, at a growth rate limited not by 
the Eddington limit but by some other physical limit on fuel 
supply which prevents accretion rates significantly in excess 
of lOM yr _1 . Such a physical limit on black-hole fuel sup- 
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ply might be imposed by accretion disc physics (e.g. the 
calculations of Burkert & Silk (2001) indicate that accretion 
disc viscosity can be expected to limit the mass consump- 
tion rate of a supcrmassivc black-hole at the centre of a 
forming spheroid to ~ 2 — 20Moyr~ 1 ) or by the physics of 
galaxy formation (e.g. Archibald et al. 2002; Granato et al. 
2001, 2003). In addition, one might speculate that the ap- 
parent lack of a substantial population of highly obscured 
supermassive black holes may be a consequence of this sub- 
Eddington accretion, as compared to the Eddington-limited 
accretion more likely to be experienced by lower mass ob- 
jects (Fabian 2003). 

6 CONCLUSIONS 

Virial black-hole mass estimates have been calculated for a 
sample of 12698 quasars in the redshift interval 0.1 < z < 2.1 
drawn from the SDSS Quasar Catalogue II (Schneider et 
al. 2003). The distribution of the quasar black- hole masses 
as a function of redshift has been presented and compared 
with the masses of dormant black-holes observed in the lo- 
cal Universe. In addition, the quasar host-galaxy properties 
implied from an application of the locally observed relation- 
ship between black-hole and bulge mass have been compared 
with the known properties of local early-type galaxies. By 
combining the black-hole mass estimates with the quasar 
bolometric luminosities the distribution of quasar accretion 
rates have been investigated. Furthermore, in combination 
with the optical quasar luminosity function, the new SDSS 
black-hole mass estimates have been used to estimate the 
number density of active black-holes at z ~ 2 as a func- 
tion of mass. Finally, the activation fraction of black holes 
at 2 ~ 2 has been estimated by comparing the estimated 
number density of active black holes with the local dormant 
black-hole mass function. The main conclusions of this study 
can be summarized as follows: 

(i) The virial black-hole mass estimates of the SDSS 
quasars are entirely consistent with an upper boundary of 
~ 3 x 10 9 Mq . This limit is consistent with both the most 
massive black-holes measured dynamically in the local Uni- 
verse, and the expected black-hole mass limit based on the 
known properties of early-type galaxies and the locally ob- 
served correlation between bulge and black-hole mass. Con- 
sequently, using the Mgll-based virial mass estimator, no 
evidence is found for a conflict between quasar black-hole 
masses in the redshift range 0.1 < z < 2.1 and the contem- 
porary, or ultimate properties of their host-galaxy popula- 
tion. 

(ii) The estimated Eddington ratios of the SDSS quasars 
show only a small level of evolution over the redshift range 
0.0 < 2 < 2.1, rising from L bo i ~ 0.15L E dd at z ~ 0.2 to 
Ltoi — 0.51/Edd at 2 ~ 2.0. Although the most luminous 
broad-line SDSS quasars are accreting at rates up to the 
Eddington limit by z ~ 1, there does not appear to be a 
significant change in accretion rates in the redshift range 
1 < 2 < 2, with the most luminous broad-line SDSS quasars 
in the current sample still accreting within a factor of ~ 2 
of the Eddington limit at z ~ 2. Consequently, it appears 
from these results that the Eddington limit is still a relevant 
physical boundary to the accretion rate of luminous broad- 
line quasars, at least at z < 2. 



(iii) Using the latest determinations of the galaxy lumi- 
nosity and stellar-velocity dispersion functions, estimates 
of the local dormant black-hole mass function for early- 
type galaxies using the Mbh ~ cr and Mbh — Ltuige re- 
lations are shown to be consistent at the high-mass end 
(Mbh > 10 s Mq). The agreement between the two meth- 
ods is found to be dependent on the adoption of similar 
levels of association scatter (~ 0.3 dex) in both correla- 
tions. Our best estimate of the total mass density of dormant 
black holes within the local early-type galaxy population is 
p bh = (2.8 ± 0.4) x 10 5 MqMpc" 3 . 

(iv) The activation fraction of supermassive black-holes 
at 2 ~ 2 is apparently an increasing function of mass, with 
an activation rate of / ~ 0.005 at M bh ~ 10 8 ' 5 Mq rising 
to / ~ 0.05 at M bh ~ 10 9 ' 5 M Q . However, it is shown that 
this result is consistent with theoretical work which predicts 
quasar lifetimes to be an increasing function of black-hole 
mass. Correcting for this expectation we find that the shape 
of the active black-hole mass function at z ~ 2 is consis- 
tent with that of the local dormant black-hole mass func- 
tion, with a normalization a factor of ~ 5 lower. Making a 
conservative correction of a factor of two to account for ge- 
ometric obscuration, the direct implication of this result is 
that the fraction of black-holes with mass > 10 8 ' 5 M which 
are in place and active at z ~ 2 is > 0.4. 

(v) A fairly robust limit on the lifetime of quasars with 
black-hole masses > 10 9 ' 5 M is found to be t Q > 2 x 10 s 
years. Black holes of this mass appear to be prevented from 
growing to masses > 10 10 Mq by some physical mechanism, 
other than the Eddington limit, which prevents accretion at 
rates > 10 Mq per year. 
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APPENDIX A: THE VIRIAL BLACK-HOLE 
MASS ESTIMATE 

The assumption that certain broad-emission lines present in 
the spectra of Type 1 quasars trace the virialized velocity 
field of the broad-line region (BLR) leads to the so-called 
virial black-hole mass estimate: 

M bh = G^RblrVblr (Al) 

where Rblr is the BLR radius and Vblr is the Keplerian 
velocity of the line-emitting gas. Recent evidence in the lit- 
erature which supports the basic assumption of virialized 
motions in the BLR has been previously discussed in intro- 
duction. In this Appendix we describe the derivation of the 
final versions of the virial mass estimators, based on both 
the H[3 and Mgll emission lines, used to estimate the SDSS 
quasar black-hole masses. Much of this material has been 
discussed previously by McLure & Jarvis (2002), and the 
reader is referred to this paper for more details. The mate- 
rial is partly repeated here for completeness, and to highlight 
a change in the adopted calibration of the Rblr — ALa re- 
lation at 3000 A which is integral to the Mgll based virial 
black-hole mass estimator. 

Al Estimating the broad-line region radius 

In the absence of reverberation mapping it is necessary to 
estimate Rblr via the correlation between Rblr and AGN 
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continuum luminosity (Wandel, Peterson & Malkan 1999; 
Kaspi et al. 2000). In McLure & Jarvis (2002) the best- 
fitting Rblr — \L\ relations were derived at both 3000 A 
and 5100 A for the 34 AGN with reverberation mapping 
estimates of Rblr from Wandel, Peterson & Malkan (1999) 
and Kaspi et al. (2000) . Using the BCES method of Akritas 
& Bershady (1996) the bisector fits at both wavelengths were 
as follows: 

Rblr = (25.2 ± 3.0) [AL 3 ooo/10 3 V] (0 ' 47±0 ' 05) (A2) 
Rblr = (26.4 ± 4.4) [AL 5100 /10 37 W/] (a61±(, ' 10) (A3) 

where Rblr is in units of light-days. The sample of 34 AGN 
with reverberation mapping estimates of Rblr is composed 
of 17 PG quasars and 17 lower-luminosity Seyfert galax- 
ies spanning the luminosity range 10 34 W< \L\ < 10 39 W. 
While this luminosity range is fairly representative of the 
SDSS quasars analysed using the Hf3 emission line (0.1 < 
z < 0.5), it is not representative of the SDSS quasars at 
z > 0.7, where we have utilized the Mgll emission line to 
estimate their black-hole masses. Simply as a result of their 
higher redshifts, these quasars reside exclusively in the lumi- 
nosity range 10 37 W< \L\ < 10 40 W. Consequently, in order 
to provide as accurate an estimate of the Rblr for each 
quasar as possible, it was decided to re-fit the Rblr — ALa 
relation at 3000 A using only those objects from the rever- 
beration mapped AGN sample with \L\ > 10 37 W. Follow- 
ing McLure & Jarvis (2002) we again adopt the BCES bi- 
sector fit as our best estimate of the Rblr. — ALa relation 
at 3000 A , with the following result: 

Rblr = (18.5 ± 6.6) [AL 3 ooo/10 37 W] (0 ' 62±0 ' 14) (A4) 

which is shown along with the previous Rblr — AZ/3000 fit 
from McLure & Jarvis (2002) in Fig Al. It can be seen 
from Fig Al that although the two fits are consistent, the 
new steeper Rblr — AL3000 relation does provide a better 
fit to the reverberation mapped PG quasars of Kaspi et al. 
with L\ > 10 37 W, which are typical of the SDSS quasars 
at z > 0.7. Consequently, it is the new steeper form of the 
Rblr — AL3000 relation which is adopted throughout the 
analysis performed in this paper, t The fit to the Rblr — 
AL5100 relation determined by McLure & Jarvis (2002) is 
adopted unchanged because of the lower luminosity range 
occupied by the SDSS quasars at z < 0.7. 

A2 Broad-line region geometry 

In Eqn Al the velocity of the BLR gas is taken as 
Vblr = fx FWHM, where FWHM is the full-width half 
maximum of your chosen broad emission line (either Hp or 
Mgll in this case) and / is a geometric factor which relates 
the FWHM to the intrinsic Keplerian velocity of the broad- 
line gas. 

Due to the fact that there is currently no consensus on 
the geometry of the BLR in radio-quiet quasars it is conven- 
tional to set / = v3/2 , which is appropriate if the orbits 
of the BLR gas are randomly orientated. Perhaps the most 

t The mean difference in black-hole mass introduced by using 
Eqn A4 in preference to Eqn A2 is only 0.086 dex and does not 
significantly change any of the results or conclusions in this paper. 




Log (AL X (3000)/W) 

Figure Al. The Rblr — AL3000 relation for the 34 AGN with re- 
verberation mapping estimates of Rblr from Kaspi et al. (2000). 
The solid line is the best-fitting BCES bi-sector fit determined for 
the full sample by McLure & Jarvis (2002); Eqn A2. The dotted 
line shows the new steeper relation adopted in this paper, which 
is a fit to those objects with XL\ > 10 37 W. 

likely alternative to randomly orientated BLR orbits is the 
possibility that the geometry of the BLR is flattened or disk- 
like in nature. In radio-loud quasars, where a rough orienta- 
tion indicator is available from the radio core-to-lobe ratio, 
there is good evidence that broad-emission line FWHM show 
a orientation dependence consistent with a disk-like geom- 
etry (e.g. Browne & Wills 1986; Brotherton 1996). Given 
the similarity between the optical emission-line spectra of 
radio-loud and radio-quiet quasars, it is perhaps not unrea- 
sonable to consider the possibility that the BLR of radio- 
quiet quasars (which dominate the SDSS sample) are also 
flattened. As discussed by McLure & Dunlop (2002), in this 
case: 

Vblr = FWHM/ (2 x sin i) (A5) 

where i is the angle between the line-of-sight and the disk 
axis (i.e. i = 0° is pole-on). It has been established for radio- 
loud quasars that the angle to the line-of-sight is consis- 
tent with being randomly distributed between 0° < i < 45° 
(Barthel 1989; Willott et al. 2000). Assuming the same sit- 
uation to be true for radio-quiet quasars, it follows that the 
mean line-of-sight angle will be < i >= 30° (averaging over 
solid angle) implying that < / >= 1. Consequently, it can 
be seen that / = 1 represents a good choice of geometric 
factor in the case of both random or flattened BLR orbits 
and is adopted throughout this paper. We note that if the 
BLR of radio-quiet quasars is flattened, objects which are 
aligned closer than 15° to the line-of-sight will have their 
black-hole masses underestimated by a factor of > 4. How- 
ever, these objects should only represent a small fraction of 
the SDSS sample (~ 10% if i is randomly distributed in the 
range 0° < i < 45°) and will likely appear as narrow-line 
objects (FWHM< 2000 km s _1 ), which are not the princi- 
pal focus of this paper. A discussion of the evidence for the 
flattened BLR hypothesis provided by the SDSS sample is 
presented in Section 3.3. 
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A3 The virial black-hole mass estimators 

Substituting the calibrations of the Rblr — ALa relations 
into Eqn Al, and making the identification that Vblr = 
Mgll or H/3 FWHM, we arrive at the virial black-hole mass 
estimators. As in McLure & Jarvis (2002), these virial esti- 
mators are then compared to the black-hole masses of the 
reverberation mapped AGN from the Kaspi et al. sample. By 
requiring that on average the new virial estimators match 
the mass estimates derived by Kaspi et al. (2000), using 
the reverberation mapping values of Rblr and rms FWHM 
measurements, small off-sets of —0.05 dex and —0.04 dex 
are introduced to the Mgll and H/3 estimators respectively. 
Accounting for these off-sets we arrive at the final versions 
of the virial black-hole mass estimators adopted throughout 
the analysis in this paper: 

M bh _ f AL30oo \°- 62 ( FWHM(MgII) \ 2 
M^-^KW^W) { krn^ ) (A6) 

M bh _ 17 f ALsioo N - 61 ( FWHM{H(3) \ 2 

M ^'{w^w) ^ kms- 1 ) 1 ' 

where the H/3 based black-hole mass estimator is unchanged 
from that derived by McLure & Jarvis (2002). 

A4 Comparing the H/3 and Mgll based estimators 

McLure & Jarvis (2002) compared the Hf3 and Mgll based 
virial estimators using a sample of 150 quasars drawn from 
the LBQS (Forster et al. 2001) and MQS samples (Baker et 
al. 1999) which had spectra covering both emission lines. 
The result of this comparison showed that Mbh(H/3) oc 
M bh (MgII) 1M)±0S)S with a la scatter in log M bh of 0.41 dex. 
Here we can re-investigate the agreement between the two 
virial mass estimators more thoroughly using the 1136 SDSS 
quasars from our sample in the redshift range 0.5 < z < 0.8 
where it was possible to estimate the black-hole mass using 
both H[3 and Mgll. The results of this comparison are il- 
lustrated in Fig A2. The top panel shows the distribution 
of the ratio of the FWHMs of the H/3 and Mgll emission 
lines. The middle panel shows the distribution of the ratio 
of broad-line region radius, as estimated via the Rblr — ALa 
correlations at 5100 A and 3000 A respectively. The distri- 
butions of the FWHM and Rblr ratios shown in the top and 
middle panels are both centred on unity, with la scatters of 
0.16 dex and 0.07 dex respectively. The bottom panel shows 
the distribution of the ratio of black- hole mass, as estimated 
via Eqn A7 and Eqn A6 respectively. The mean log ratio of 
the black-hole mass estimates is 0.013 with a la scatter of 
0.33 dex. In conclusion, this comparison demonstrates that 
the Mgll and H/3 based virial mass estimators return the 
same mass estimates within a factor of ~ 2 on average, and 
that statistically are entirely interchangeable. This conclu- 
sion is re-enforced by the smooth transition between the two 
estimators at z ~ 0.7 shown in Fig 1. 



APPENDIX B: MODEL FITTING 

Due to the size of the sample of SDSS quasars analysed 
here, the fitting of the H/3 and Mgll emission lines and the 
determination of the continuum luminosities at 5100 A and 
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Figure A2. The top panel shows a histogram of the log of the 
ratio of H/3/MgII FWHM for 1136 objects in the redshift interval 
0.5 < z < 0.8. The mean log ratio is -0.004 ± 0.005, with a 
dispersion of 0.16 dex. The middle panel shows a histogram of 
the log of the ratio of Rblr as estimated via the Rblr — ALa 
relation at 5100 A and 3000 A respectively. The mean log Rblr 
ratio is 0.0010±0.0021, with a dispersion of 0.07 dex. The bottom 
panel shows a histogram of the log of the ratio of virial black-hole 
mass estimate based on the H/3 and Mgll emission lines and the 
Rblr — <^A relations at 5100A and 3000A respectively. The 
mean log ratio is 0.013 ± 0.010 with a dispersion of 0.33 dex. 
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Figure Bl. An example of the emission-line fitting process in the HP region for an SDSS quasar at z=0.17. The upper plot shows 
the original spectrum shifted to the rest-frame (thin line) with the best-fitting model over-plotted as the thick line. Also shown in the 
upper plot is the fitted power-law quasar continuum (dotted line). The lower plot shows the original spectrum after subtraction of the 
best-fitting quasar continuum and iron-emission template. The thick solid line shows the best-fitting model, while the thin solid line 
shows the individual emission-line components which comprise the best-fitting model. 
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Figure B2. An example of the emission-line fitting process in the Mgll region for an SDSS quasar at z=0.88. The upper plot shows 
the original spectrum shifted to the rest-frame (thin line) with the best-fitting model over-plotted as the thick line. Also shown in the 
upper plot is the fitted power-law quasar continuum (dotted line). The lower plot shows the original spectrum after subtraction of the 
best-fitting quasar continuum and iron-emission template. The thick solid line shows the best-fitting model, while the thin solid line 
shows the individual emission-line components which comprise the best-fitting model. 



3000 A was undertaken in a fully automated fashion. The 
various stages in the emission-line modelling of both Hfl 
and Mgll are discussed below. 



Bl Continuum and iron-template fitting 

In order to reliably measure the line-widths of the H/3 and 
Mgll emission lines it is necessary to subtract the surround- 
ing iron emission. This subtraction process can often be im- 
portant for fitting the Hj3 emission line, but is crucial for 
fitting the Mgll emission line where the wings are usually 
strongly blended with the surrounding Fell emission. The 



approach adopted here for dealing with the iron emission is 
a variant on the technique developed by Boroson & Green 
(1992), which is based on the fitting of an iron-emission 
template derived from the spectrum of narrow-line Seyfert 
galaxy (NLS1) Izwl. As a member of the NLS1 class, Izwl 
is ideal for this purpose because it displays extremely strong 
iron emission around both H/3 and Mgll with a FWHM of 
~ 900 km s _1 . Consequently, the strength and FWHM of an 
iron template based on Izwl can be adjusted to match the 
iron emission of broad-line quasars where the FWHM of the 
iron emission is > 2000 kms" 1 . The iron emission templates 
used in the modelling are taken from two different sources, 
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but are both derived from spectra of Izwf. For the mod- 
elling of the Hf3 emission line we have used the optical Fell 
emission template derived by Marziani et al. (2003), which 
is made freely available by the authors. For the modelling of 
the Mgll emission line we have reconstructed the Fell+Felll 
emission template derived by Vestergaard & Wilkes (200f ) 
from the data published within that paper. 

The first stage in the modelling process of both the 
HP and Mgll emission lines is the simultaneous fitting of 
the underlying quasar continuum and the iron-emission tem- 
plate, where the quasar continuum is assumed to be a simple 
power- law (/„ oc v~ a ). For both emission lines this fitting 
was performed in two wavelength windows, which conserva- 
tively excluded the Hp and Mgll emission regions to avoid 
biasing the iron-emission fit with the wings of the line emis- 
sion. For the H/3 emission line the continuum+Fell fitting 
windows were 4435 A -> 4700 A and 5100 A ->5535 A . The 
equivalent regions for the Mgll emission line where 2f 00 A — > 
2675 A and 2925 A ^3090 A . During the fitting process the 
normalization and slope of the continuum were left as free 
parameters, as were the luminosity and FWHM of the iron 
template. Upon completion of this process the best-fitting 
continuum+iron combination was subtracted from the spec- 
tra, leaving the isolated Mgll or _ff/3+[Olll] emission lines. 
The fitting of the Mgll and H/3 emission lines then pro- 
ceeded in slightly different fashions, and are therefore de- 
scribed separately below. 

B2 Fitting the H/3-\-[Olll] emission-line complex 

The H/3 and [Oni] emission lines were modelled simultane- 
ously as a combination of four components: broad and nar- 
row H/3 plus narrow components for [Oni] 4959 A and [Oni] 
5007 A . In order to limit the degrees of freedom involved in 
the fitting several restrictions were invoked. Firstly, the cen- 
tral wavelengths of the narrow Hp component and the [Oni] 
4959 A component were constrained to lie at their labora- 
tory wavelength offsets with respect to [Oni] 5007 A. This 
meant that there were only two free central wavelengths 
involved in the fit, that of the broad HP component and 
[Om] 5007 A . Secondly, the FWHM of the narrow Hp com- 
ponent and the two [Oni] components were constrained to 
have the same value. The final restriction imposed was that 
the FWHM of the narrow components was constrained to 
be less than 2000 km s _1 . 

Four different representations of the line profiles were 
tried for each object. Firstly, two fits were performed where 
all the line profiles (broad+narrow HP and [Om] 4959 
& 5007) were first modelled as gaussians and then as 
lorentzians. Secondly, two further fits were performed where 
the broad Hp profile was first modelled as a gaussian, with 
all the narrow profiles modelled as lorentzians, and then sec- 
ondly with broad HP modelled as a lorentzian with all the 
narrow profiles modelled as gaussians. The FWHM value 
used in the virial mass estimates is that of the best-fitting 
model profile, in terms of \ 2 1 within the line-fitting region 
(4700 A -> 5100 A). 

B3 Fitting the Mgll emission line 

The first stage of the modelling of the Mgll emission line 
involved the fitting of a two-component, broad+narrow pro- 



file. This proceeded in similar fashion to the modelling of 
the HP emission line, with both the narrow and broad com- 
ponents first represented as either gaussians or lorentzians. 
In addition, two further model fits were tried for each ob- 
ject, the first comprising a gaussian broad component and 
a lorentzian narrow component, and the second a lorentzian 
broad component and a gaussian narrow component. During 
this stage of the modelling the central wavelength, FWHM 
and flux of both the narrow and broad components were 
treated as free parameters. The only restriction invoked was 
that the narrow component had a FWHM < 2000 km s _1 . 
Both the narrow and broad components were modelled as 
doublets separated by their laboratory wavelength differ- 
ence, and all of the emission-line models were fitted within 
the wavelength range 2700 A -» 2900 A . 

Unlike the HP emission line, which often displays a clear 
point of inflection between the narrow and broad compo- 
nents, with Mgll it is often unclear whether a significant, 
separate narrow component is present. To account for this, 
the Mgll emission line of each quasar was also modelled as 
a single broad lorentzian or gaussian doublet. For the pur- 
poses of the virial mass estimates the FWHM of the broad 
component of the best-fitting two-component model profile 
was adopted as the best FWHM estimate for each quasar, 
provided that two criteria were satisfied. Firstly, the A\ 2 be- 
tween the best-fitting two-component and best-fitting single- 
component models was required to be significant at the 3a- 
level, i.e. the extra degrees of freedom in the two-component 
model had to be justified statistically. Secondly, the equiv- 
alent width (EW) of the broad component was required to 
be greater than two-thirds of the total Mgll EW, ensuring 
that the FWHM of the broad component was only adopted 
if it constituted the dominant fraction of the line flux. If 
either of these two criteria were not satisfied the FWHM of 
the best-fitting single-component model profile was adopted. 
The typical uncertainty in the fitted FWHM values for both 
the H/3 and Mgll emission lines is ±15%. 



APPENDIX C: LUMINOSITY CALIBRATION 

For the purposes of the analysis carried out within the 
paper it was necessary to provide a consistent estimate of 
the quasar bolometric luminosities, based on the continuum 
close to either the Hp or Mgll emission lines. Furthermore, 
in order to perform the estimate of the quasar activation 
rates in Section 5, it was necessary to provide a calibra- 
tion between the quasar bolometric luminosities estimated 
via the spectra modelling, and the absolute B— band mag- 
nitudes used in the Boyle et al. (2000) determination of the 
2dF+LBQS quasar luminosity function. Both of these cali- 
bration issues are discussed below. 

CI Estimating the bolometric luminosity 

The bolometric luminosities of each of the SDSS quasars are 
estimated from their monochromatic \L\ luminosity at ei- 
ther 5100 A or 3000 A , depending on their redshift. The ap- 
propriate bolometric correction factors have been estimated 
by extrapolating from the bolometric correction at 2500 A 
determined by Elvis et al. (1994), employing their median 
quasar SED. Taking into account the recent re-evaluation of 
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which is shown in Fig CI. The la scatter at a fixed abso- 
lute magnitude is 0.14 dex, as highlighted by the dashed 
lines. Three of the 372 objects were found to be large 
outliers and were excluded from the fitting process. The 
most extreme outlier (shown as the filled circle in Fig CI) 
is the z = 1.296 quasar SDSS J133855.04-010229.9 (2QZ 
J133855. 0-010230). In this case, the source of the discrep- 
ancy appears to be an erroneous Bj magnitude listing in 
the 2dF 10K catalogue. The Bj magnitude for this quasar 
is listed as Bj — 20.82. However, the SDSS photometry 
for this object finds a g— band magnitude of 18.42 with a 
g — r colour of 0.18. Even accounting for the switch between 
AB/Vega magnitudes and galactic extinction, it is clear that 
the SDSS and 2dF photometry for this object are inconsis- 
tent, whereas the SDSS PSF magnitude and flux calibrated 
spectra are consistent. 



Log (WW) 



Figure CI. The calibration (solid line) between absolute 
B— band magnitude and estimated bolometric luminosity, based 
on N=372 quasars common to both the SDSS Quasar Catalog II 
and the 2dF 10K quasar catalogs. The two dashed lines delineate 
the lcr = 0.14 scatter around the best-fit. The outlying object 
shown as a filled circle is suspected of having an incorrect mag- 
nitude in the 2dF 10K quasar catalogue (see text for discussion). 

the average 2500 A^ 2 keV X -ray slope (Elvis, Risaliti & 
Zamorani 2002) we estimate the bolometric correction fac- 
tors to be 5.9 and 11.3 at 3000 A and 5100 A respectively. 

Due to the fact that the SDSS sample of quasars anal- 
ysed here includes 1136 objects in the redshift range 0.5 < 
z < 0.8, whose spectra cover both 3000 A and 5100 A, we 
have an opportunity to verify that our choice of bolometric 
correction factors lead to consistent bolometric luminosity 
estimates at both wavelengths. Using correction factors of 
5.9 and 11.3 we find that the mean log ratio of bolometric 
luminosity estimates is : < log(i(, o ;(5100)/I/6 O ;(3000)) >= 
0.062 ± 0.003. This result is inconsistent with the expected 
ratio of unity, and is equivalent to the 5100 A based bolo- 
metric luminosity estimate being 15% greater on average. 
Consequently, we chose to revise downward the 5100 A bolo- 
metric correction factor to 9.8, which ensures a mean ratio 
of unity for the 1136 SDSS quasars whose spectra cover both 
5100 A and 3000 A , with a la scatter of less than 30%. 



C2 The Mb - bolometric luminosity calibration 

The calibration between absolute B— band magnitude and 
the bolometric luminosities derived from the fits to the flux 
calibrated SDSS spectra was done using a sample of 372 ob- 
jects common to the 2dF 10K quasar catalogue (Croom et 
al. 2001) and the new SDSS Quasar Catalogue II (Schnei- 
der et al. 2003). The absolute B— band magnitudes for each 
quasar have been calculated from the apparent Bj magni- 
tudes listed in the 2dF 10K quasar catalog, using the K- 
corrections and quasar colours from Cristiani & Vio (1990) 
for consistency with Boyle et al. (2000). The best-fitting re- 
lation to the correlation between bolometric luminosity and 
absolute B— band magnitude is: 



M B = -2.66(±0.05) log [L bol /W] + 79.36(±1.98) (CI) 
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